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Abstract 

We consider settings in which T multi-antenna transmitters and K single-antenna receivers concur- 
rently utilize the available communication resources. Each transmitter sends useful information only to 
its intended receivers and can degrade the performance of unintended systems. Here, we assume the 
performance measures associated with each receiver are monotonic with the received power gains. In 
general, the systems' joint operation is desired to be Pareto optimal. However, designing Pareto optimal 
resource allocation schemes is known to be difficult. In order to reduce the complexity of achieving 
efficient operating points, we show that it is sufficient to consider rank-1 transmit co variance matrices 
and propose a framework for determining the efficient beamforming vectors. These beamforming vectors 
are thereby also parameterized by T(K — 1) real-valued parameters each between zero and one. The 
framework is based on analyzing each transmitter's power gain-region which is composed of all jointly 
achievable power gains at the receivers. The efficient beamforming vectors are on a specific boundary 
section of the power gain-region, and in certain scenarios it is shown that it is necessary to perform 
additional power allocation on the beamforming vectors. Two examples which include broadcast and 
multicast data as well as a cognitive radio application scenario illustrate the results. 

Part of this work has been performed in the framework of the European research project SAPHYRE, which is partly funded 
by the European Union under its FP7 ICT Objective 1.1 - The Network of the Future. This work is also supported in part by 
the Deutsche Forschungsgemeinschaft (DFG) under grant Jo 801/4-1. 

The authors are with the Department of Electrical Engineering and Information Technology, Dresden University of Technology, 
01062 Dresden, Germany. E-mail: {Rami.Mochaourab,Eduard.Jorswieck}@tu-dresden.de. Phone: +49-351-46332239. Fax: +49- 
351-46337236. 

Part of this work has been presented at IEEE International Zurich Seminar on Communications, Zurich, Switzerland, March 
3-5, 2010 |1| and International Workshop on Cognitive Information Processing, Elba Island, Italy, June 14-17, 2010 (2). 



November 2, 2010 



DRAFT 



2 



I. Introduction 

Interference is known to be one of the major factors that limits the performance of a communication 
system in a wireless network. This situation is common in multiuser settings when the systems concur- 
rently share the available communication resources. In general interference networks, the performance 
measure of individual users is described by a utility function. This function depends in a monotonic 
way on the received signal power, interference signal power and noise power. The joint operation of the 
systems is efficient if it is not possible to improve the performance of one system without degrading 
the performance of another. In this case, the operating point is said to be Pareto optimal. It is always 
desired to design resource allocation schemes that lead to Pareto optimal operation points. In this way, 
the available communication resources are utilized efficiently to grant efficient operation of the systems. 
However, developing efficient resource allocation schemes is not straightforward and proves to be difficult. 
For instance, the problem of finding the maximum sum-rate or the proportional-fair operating point in 
the multiple-input single-output (MISO) interference channel (IC) is proven to be strongly NP-harcQ (3). 

In a multiuser setting, efficient operation of the systems requires the transmitters to maximize the 
power gain at intended receivers and also minimize the power gain at unintended receivers. In this work, 
we characterize the transmission strategies of each transmitter that are relevant to achieve Pareto optimal 
operating points. Moreover, we parameterize these by real values between zero and one. In this way, the 
set of efficient transmission strategies is confined and represented by low dimensional real parameters. 
This result tremendously reduces the complexity of designing efficient resource allocation schemes, and 
the parametrization can be utilized for low complexity coordination between transmitters. 

We give a brief reference to related work in the MISO IC and a few of their applications. The MISO IC 
is an example of an interference network where the systems consist of transmitter-receiver pairs. For the 
two-user case, real-valued parametrization of each transmitter's efficient beamforming vectors is provided 
in H. The beamforming vectors that achieve Pareto optimal points are proven to be a linear combination 
of zero-forcing (ZF) transmission and maximum ratio transmission (MRT). Based on this characterization, 
a monotonic optimization framework is developed in @ to find maximum sum-rate, proportional-fair 
and minimax operating points. The parametrization in [fH relates to a parametrization using the virtual 
SINR framework in Q. The use of this framework is motivated by the design of distributed algorithms 
that require local channel state information (CSI) at each transmitter. This framework is extended to 
the precoding design in MIMO settings in Q. The concept of combining the MRT and ZF strategies 

'interestingly, these problems are efficiently solvable if rate requirements or interference constraints on each system are fixed. 
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is important for developing so-called distributed bargaining algorithms. These algorithms improve the 
operation of the systems from the noncooperative outcome (H, 0. In (H, a distributed bargaining 
algorithm is developed which requires one bit signaling between the transmitters. Extension to the 
precoding design in the multiple-input multiple-output (MIMO) case is given in [10]. In [9], a similar 
distributed beamforming algorithm in the MISO IC is proposed for the case of statistical CSI at the 
transmitters. Also utilizing the parametrization in |4), a distributed bargaining process is proposed in 
iTTTi which requires four bit signaling between the transmitters. The process is proven to converge to an 
operating point arbitrarily close to the Pareto boundary and dominates the noncooperative outcome of the 
systems. In lfl2l . the high signal to interference plus noise ratio (SINR) approximation of the achievable 
sum-rate of a system pair is utilized to determine suboptimal joint transmission strategies. The achieved 
performance is shown to be better than the joint MRT and joint ZF strategies. 

In the K-user MISO IC, complex-valued parametrization of the Pareto boundary of the MISO IC 
rate -region is derived in |[T3l which requires K(K — 1) complex-valued parameters in order to attain all 
Pareto optimal points. In [14], the i^-user MISO IC is considered with the capabilities of time sharing 
the resources between the links. All points on the Pareto boundary of the MISO IC rate region are 
achieved with K(K — 1) real valued parameters each between and n. In lfl5l . the authors characterize 
the Pareto boundary of the MISO IC through controlling interference temperature constraints (ITC) at 
the receivers. Each Pareto optimal rate tuple is achieved iteratively when each transmitter optimizes its 
transmission constrained by the ITCs. It is shown that K(K — 1) real valued parameters, each between 
zero and a value depending on the channel vectors, are needed to achieve all Pareto optimal points. ITC 
is a terminology used in cognitive radio scenarios under the underlay paradigm [16]. It quantifies the 
amount of interference from the secondary transmitters that is tolerated by the primary users. 

In ifTTI , joint linear precoding is investigated taking into account the signaling overhead between the 
transmitters. The rate-region achieved with joint precoding is larger than the MISO IC rate -region, and 
all Pareto optimal beamforming vectors are parameterized by K(K — 1) complex-valued parameters. For 
the same setting, a recent result in [18 ] reduces the number of parameters to K + L real- valued scalars, 
each between zero and one, where L is the number of linear constraints on the transmission. Linear 
precoding MIMO IC algorithms are moreover investigated in |fl9l for a two-user system. 

While the above mentioned results are provided for the MISO IC setting, the parametrization of efficient 
transmission strategies in a general multiuser setting is not straightforward. Moreover, neither the ITC- 
based |fT5l nor the Lagrangian-based |[T4l . |fT8l characterizations can be generalized to our framework. 
A further example of a MISO multiuser setting which can be applied to our framework is when a 
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single transmitter sends common information to K single-antenna receivers. This setting corresponds 
to multicast transmission. Since the transmission rate depends on the weakest link in the system, the 
transmitter optimizes its transmission to achieve max-min-fairness at the receivers |20l . The multicast 
beamforming problem to achieve max-min-fairness is proven to be NP hard for K > N 11211 . where N 
is the number of transmit antennas. In |22], the two-user multicast max-min-fair problem is studied, and 
the set of beamforming vectors which includes the solution of the max-min-fair problem is characterized. 

In this work, we provide a general framework for parameterizing the transmission strategies of each 
transmitter which are relevant to achieve Pareto optimal points. This framework is applicable to settings 
where the utility functions of the systems are monotonic in the received power gains. The contributions 
and outline of this paper are as follows: 

• We investigate the properties of efficient transmission of a single transmitter. These properties are 
acquired on studying the transmitter's power gain-region (Section ITTT1 ). The power gain-region is 
composed of all jointly achievable power gains at the receivers. Of interest are the transmission 
strategies which achieve its boundary part in a specific direction. We prove that the boundary of 
the power gain-region is convex and always achieved with single-stream beamforming (Lemma [3]). 
Due to these properties, the corresponding strategies are characterized by real-valued parameters 
(Theorem Q]). Furthermore, we characterize under which conditions power control is needed for 
efficient transmission, (i) When the number of transmit antennas is greater than or equal to the 
number of receivers K (Section [III- At . we prove that full power transmission achieves all boundary 
points (Lemma[2]). In this case, K — 1 real- valued parameters, each between zero and one, are needed 
to parameterize the beamforming vectors, (ii) When the number of transmit antennas is strictly less 
than the number of receivers (Section IIII-BI ). we characterize the transmission strategies for which 
power control is needed. For this case, an additional real-valued parameter between zero and one is 
needed that varies the power level at the transmitter. 

• We utilize the developed single-transmitter framework for the multiple-transmitter case (Section [TV). 
Based on the network setting and the monotonicity properties of each receiver's utility function, the 
boundary part which is relevant for Pareto optimal operation is determined for each transmitter's 
gain-region. Consequently, each transmitter's efficient strategies are parameterized (Theorem [2]). We 
provide an example setting which includes broadcast and multicast data, and we apply the developed 
framework to this setting (Section II V- Al l. Moreover, we apply the framework to the A'-user MISO 
IC (Section [V- Al l. As a special case, the result for the two-user MISO IC in |fT3l follows. In addition, 
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we give an alternative characterization of the efficient transmission strategies (Corollary [T]) which 
is motivated by the application of null-shaping constraints in underlay cognitive radio scenarios 
(Section IV-BI ). We prove that all Pareto optimal operating points can be characterized through the 
design of null-shaping constraints on noncooperative secondary transmitters. Extensions to the case 
of multiple antennas at the receivers is covered in Section [VT] 
Notations: Column vectors and matrices are given in lowercase and uppercase boldface letters, re- 
spectively. The notation Xk,e describes the £th component of vector xy-. The Euclidean norm of a vector 
a, a € C j is written as ||a||, and the absolute value of b, b G C, is \b\. (-) H denotes the Hermitian 
transpose. The ith eigenvalue of a matrix Z is denoted by \1%{Z). The eigenvector which belongs to the 
ith eigenvalue of the matrix Z is denoted by vi(Z). We always assume that the eigenvalues are ordered in 
nondecreasing order such that ^i{Z) < Hi + \(Z). Moreover, the eigenvectors corresponding to the largest 
and smallest eigenvalues of a matrix Z are specified as v m;ix (Z) and v m i n (Z), respectively. The notation 
Z y means that Z is positive semidefinite. The rank and trace of a matrix Z are given by rank (Z) and 
tr(Z), respectively. The orthogonal projector onto the column space of Z is Ylz ■= Z(Z H Z) 1 Z H . 
The orthogonal projector onto the orthogonal complement of the column space of Z is := I — Hz, 
where i" is an identity matrix. E(-) denotes statistical expectation. The set of non-negative real numbers 
is R + . The cardinality of a set K is written as |/C|. 



We consider T transmitters and K receivers sharing the same spectral band. Define the set of trans- 
mitters as T := {1, T} and receivers as K := {1, K}. Each transmitter sends useful information to 
at least one receiver. For transmitter k,k G T, let K{k) C JC denote the set of its intended receivers for 
which useful information is sent to, and let JC(k) = K\K,{k) be the set of its unintended receivers. Each 
transmitter k is equipped with N k antennas, and each receiver with a single antenna. The quasi-static 
block flat-fading instantaneous channel vector from transmitter k, k G T, to receiver i, I G /C, is denoted 
by h^i G C fcXl . The transmit covariance matrix of transmitter k is given as Q k G <C NkXNk , Q k >z 0. 
We do not make any assumptions on the number of data streams applied at the transmitters. The basic 
model for the matched-filtered, symbol-sampled complex baseband data received at receiver £ is 



where is the symbols vector transmitted by transmitter k and ri£ are the noise terms which we model 
as independent and identically distributed (i.i.d.) complex Gaussian with zero mean and variance a 2 . 



II. System and Channel Model 



T 




(1) 



k=l 
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Each transmitter has a total power constraint of P := 1 which leads to the constraint tr (Q k ) < 1, A £ T. 
Throughout, we define the signal to noise ratio (SNR) as I /a 2 . The feasible set of co variance matrices 
for transmitter k is defined as 

S k := {Q k e C N » xN * : Q k y 0,tr(Q fc ) < l}. (2) 

Note that S k is compact and convex. We assume each transmitter has local CSI, i.e., it has perfect 
knowledge of the channel vectors only between itself and all receivers Q. This ideal scenario serves as 
an upper bound to the more realistic case in which imperfect or partial CSI at the transmitters is available. 
Extensions in this direction are reported in ||23"1 . |24|. In these works, Pareto efficient transmission 
strategies are characterized for the two-user MISO IC with partial CSI at the transmitters. 

A. Assumptions on Performance Measure 

The performance measure of a system in an interference network is usually described by a utility 
function. The utility function associated with a receiver depends on the power gains originating from the 
transmitters in the network. Define the power gain achieved by transmitter k at a receiver I as 

Xk,e(Qk) = h ke Qkh ke , (3) 

where xi{Q k ) € M+ since Q k is positive semidefinite. The utility function associated with a receiver i 
is defined as ug : — >• R + , where T is the number of transmitters in the network. 
Assumption 1: The utility function ug,£ G fC, has the following properties: 

A. If £ £ fC(k), then U£ is monotonically increasing in the power gain from transmitter k, i.e., 

^Ol,£(Ql); -,X T ,£(Qt)) < u i{%w{Ql), -iX k ,i{Q k ), ...,XT,e(Q T )), (4) 

for x kjl {Q x ) < x k/ (Q k ). 

B. If £ £ K{k), then ue is monotonically decreasing in the power gain from transmitter k, i.e., 

ue( x iAQi)' ■■■i x t,i(Qt)) > uifxxjiQx), ...,x k/ (Q k ), ...,x T> i(Q T )\ (5) 

for x k: i(Q k ) < x k ,i{Q k ). □ 
Assumption [U describes the settings where the performance measure at a receiver increases monoton- 
ically with increased power gain from intended transmitters and decreases monotonically with increased 
power gain from unintended transmitters. An example utility function which satisfies Assumption Q] is 
the signal to interference plus noise ratio (SINR). 
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The utility region is the set of all achievable utility tuples defined as: 



U := {(iti(sci,i(Qi), • • .,x T ,i(Q T )), ■ ■ ■ ,uk(xi,k(Qi), ■ ■ -,%T,K (Q T )))-.Q k eS k ,keT}cRf. (6) 



The efficient operating points in the utility region correspond to those in which it is impossible to improve 
the performance of one system without simultaneously degrading the performance of at least one other 
system. Such operating points are called Pareto optimal and are defined formally as follows. 

Definition 1: A tuple (ui, ...,uk) G U is Pareto optimal if there is no other tuple {u'\, ...,u'k) G U 
such that (it'i, ...,u'k) > (ui, ...,uk), where the inequality is component-wise and strict for at least one 
component. The set of all Pareto optimal operating points constitutes the Pareto boundary {VB) of U. □ 

Next, we give an example setting where the systems' utility functions satisfy Assumption [TJ 
B. Example Setting 

Consider two transmitters each using three transmit antennas, and three single antenna receivers as 
depicted in Fig. Q] The operation of the systems is as follows: 

• Broadcast Channel (BC): Transmitter 1 transmits different useful data to receivers 1 and 2 simulta- 
neously. We assume transmitter 1 chooses the transmit covariance matrices Q u with tr (Qn) = Pn 
for receiver 1 and Q 12 with tr (Q 12 ) = P12 f° r receiver 2. Hence, transmitter 1 can be considered 
as two virtual transmitter^, 11 and 12, coupled by the total power constraint, pn + pu < 1. The 
receivers are identified in the following receiver sets: 1 G /C(ll), 1 € /C(12), 2 G /C(12), 2 G /C(ll). 

• Multiple Access Channel (MAC): Transmitters 12 and 2 send distinct useful information to receiver 
2. Receiver 2 decodes the data from transmitter 12 and 2 successively. Thus, 2 G /C(12),2 G JC(2). 

• Multicast: Transmitter 2 sends common useful data in a multicast to receivers 2 and 3. The receivers 
are identified in the following receiver sets: 2 G /C(2), 3 G /C(2). 

• Interference Channel (IC): Transmitter 2 induces interference on receiver 1, while transmitter 1 
induces interference on receiver 3. 

The receiver sets are summarized in Fig. [T] and the solid and dashed arrows refer to useful and not 
useful signal directions, respectively. The achievable rate at receiver 1 is 



ui(x 11j1 (Q 11 ),x 1 2,i(Q 1 2),X2,i(Q 2 )) = log 2 1 + 



a 



2 + hf 1 Q 12 h n + hf 1 Q 2 h 2 i 



h u Q u hii 




2 This transmission strategy is suboptimal, however less complex and more robust than dirty paper coding 1251 . 
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Fig. 1. An example setting for the described system model. There exist two transmitters, each equipped with three antennas, 
and three single antenna receivers. The solid arrows refer to the intended receivers of a transmitter, while the dashed arrows 
refer to interference directions. 



which is monotonically increasing in xn i i(Q 11 ) and monotonically decreasing in the power gains from 
transmitters 12 and 2. The utility at receiver 2 is its sum capacity |[26l . 

I lr\ \ m \ \\ 1 fi , ^VlQ 12^12 + h 22 Q 2 h 22 \ 
«2(a:il,2(Qll), Zl2,2(Ql 2 ), ^2, 2 (Q 2 )) = lo S2 1 + 7, , H ^ , , (8) 

V cr 2 + hf 2 Q u h 12 J 

which is monotonically increasing in 0:12,2(^12) an ^ ^2,2(^2)- The utility function at receiver 3 is the 
achievable ratd_, 

u 3 (x n>3 (Q n ),x 12i3 (Q 12 ),x 2!3 (Q 2 ))=log 2 (l+ u h^u 2H2 \h„ u )■ (9) 

V a 1 + hf 3 Q u h 13 + hf 3 Q l2 h 13 J 

which is monotonically increasing in x 2:3 (Q 2 ). The utility functions in ©-(O satisfy properties A and B 
in Assumption Q] We return to this setting in Section IIV-AI after we formalize the framework for efficient 
beamforming and resource allocation. 

III. Power Gain Region 

In this section, a single transmitter k,k G T, is considered along with all K receivers. The subscript 
k in all terms referring to the single transmitter is omitted for convenience. For example, hi is written 
instead of h^. The receiver sets of the transmitter are written as K, and K instead of JC{k) and JC(k). 

3 Note that the transmission rate at transmitter 2 has to be chosen such that both receiver 2 and 3 can decode the data 
successfully. We do not consider this requirement in {8} and l|9} since this is beyond the scope of this paper. However, these 
rates can be achieved using rateless coding |27|, |28|. 
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We return to include the indication to the transmitter in the next sections when multiple transmitters 
are considered. Here, we study the transmission effects of a single transmitter on all existing receivers. 
Thereby, we characterize its transmit covariance matrices that are relevant for its efficient operation in the 
multiuser system. For all feasible transmit covariance matrices, a power gain-region of a single transmitter 
is the set that consists of all joint power gains achievable at the receivers. The power gain-region of a 
transmitter is defined as 

n := {(ari(Q), ...,x K (Q)) : Q G S} C , (10) 

where S is defined in (fSJ). Note that x\(Q), ...,xk{Q) are the main elements of the utility functions in 
©-(O. An important property of the gain-region fi is its convexity. Having this property is convenient for 
characterizing the points on its boundary using simple programming problems based on the Hyperplane 
Separation theorem |29l Theorem 1.3]. 

Lemma 1: The set Q in ( fTOl ) is a compact and convex set. 

Proof: The proof is provided in Appendix |Al ■ 

Points that lie on the boundary of the gain-region Q are of interest since these points characterize 
extreme power gains achievable at the receivers. At these points, the transmitter cannot increase the 
power gain in one direction of the gain-region without decreasing the power gain in any other direction. 
We give an example later in this section which further clarifies the importance of the boundary points 
of the gain-region for efficient operation of a transmitter. Next, we formalize the boundary of the set 
following the definitions in ||5]. There, these definitions were used to derive the solution of monotonic 
optimization problems [30], |29l . 

Definition 2: A point y G Wl is called upper boundary point of a compact convex set C if y G C 
while the set |y' G R+ : y' > y] C M™ \ C, where the inequality in y' > y is componentwise. The set 
of upper boundary points of C is denoted by d + C. □ 

Definition |2] describes only one boundary part of a compact convex set. The straightforward extension 
to describe all boundary parts of this set is to define its upper boundary in direction e,e G {— 1, +l} n . 
For this purpose, we first need the following definition. 

Definition 3: A vector x dominates a vector y in direction e, written as x > e y, if xgeg > ygeg for 
all I, 1 < I < n, and the inequality has at least one strict inequality. □ 

Definition 4: A point y G R" is called upper boundary point of a compact convex set C in direction 
e if y G C while the set {y' G R™ : y' > e y}cK™\ C. We denote the set of upper boundary points in 
direction e as d e C. □ 
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Fig. 2. An illustration of a two-dimensional gain-region and its upper boundaries in directions ei = [1, 1], e.i = [1,-1], and 
es = [-l,l]. 

An illustration of a two-dimensional power gain-region is given in Fig. [2] The direction vectors e\, e2, 
and e3 refer to three different parts of the boundary. For the choice e = 1, the upper boundary in direction 
e is the usual upper boundary as in Definition [2j i.e., d + Q, = d 1 ^}. 

We draw an example to illustrate the importance of the boundary parts of the gain-region. Assume as 
in Fig. |2j there exist two single antenna receivers and a transmitter. Assume, receiver 1 is the intended 
receiver of the transmitter, i.e. 1 = K, and receiver 2 is its unintended receiver such that 2 = K. For 
efficient operation in the setting, the transmitter is interested in maximizing its power gain at receiver 
1 and also interested in minimizing the power gain on receiver 2. If we seek to characterize the set of 
efficient transmission strategies of the transmitter, then the boundary part corresponding to e2 = [1,-1] is 
relevant. Any power gain tuple which is inside the power gain-region and does not lie on its boundary is 
not relevant, since the transmitter can increase the power gain to its intended receiver without changing 
the power gain to the unintended receiver. The direction vector e.2 corresponds to the boundary part 
which includes the maximum achievable power gain at receiver 1 and also the minimum achievable 
power gain at receiver 2. These extreme points correspond to MRT and ZF transmission strategies, 
respectively. According to Fig. [2l MRT achieves ||/ii|| 2 power gain at the first receiver and x\ at the 
second receiver. On the other hand, ZF transmission achieves zero power gain at the second receiver and 
x\ at the first receiver. From this example, it can be observed that the boundary part which corresponds 
to the transmitters efficient strategy set is the one where the direction vector has positive component 
corresponding to the intended receiver and negative component corresponding to the unintended receiver. 
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Since a single transmitter has at least one intended receiver, otherwise it will not operate, the direction 
vector with all components equal to — 1 is not of interest. Therefore, we define the set of feasible directions 



Next, we study the transmission strategies that achieve the boundary points in c? e J7, e 6 £. 

The relation between the number of existing receivers and the number of available antennas at the 
transmitter is important to distinguish whether power control is needed to achieve all boundary parts of 
the power gain-region. First, we consider the case where the number of transmit antennas is greater than 
or equal to the number of existing receivers K. For this case, a comprehensive study of the gain-region 
is given due to which we gain a link and some insights to a mathematical field of research in matrix 
analysis. Afterwards, we consider the case in which the number of antennas is strictly less than K. This 
case is addressed more briefly since the tools needed for the analysis are similar to those in the first case. 

A. The number of transmit antennas satisfies N > K 

In this section, we assume N > K and the channel vectors to the receivers are independently distributed. 
These assumptions lead to linear independence of the channel vectors with probability one. In this case, it 
is possible to achieve power gain on one receiver and simultaneously null the power gain at the remaining 
receivers. Therefore, the gain-region has boundary points that lie on each axis as the points x\ and xJ, m 
Fig. |2l As a result, all upper boundary points of the gain-region are achieved with full power. In case the 
number of antennas at the transmitter is strictly larger than K, the gain-region has part of its boundary 
the points between x\ and the origin. 

Lemma 2: Transmit covariance matrices from the set 



Lemma |2] supports the result in |[3ll . where it is shown that all Pare to efficient operating points in the 
MISO IC correspond to full power transmission when the number of antennas is larger than or equal to 
the number of receivers. The power gain-region achieved with full power transmit covariance matrices 
is defined as 



£:={-l,ir\{-l} 



(11) 



S :={Q :Qh 0,tr(Q) = 1} 



(12) 



achieve all points in d e Q, e € £. 

Proof: The proof is provided in Appendix |B] 




(13) 
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Since the power gain at a receiver £ can equivalently be formulated as xg(Q) = hfQhi = tr [Qh(hf\, 
the set Q, in (fT3l ) is rewritten as 



This set is referred to in 
$7 is compact and convex 



n = {(tr(Qfrifrf) ,...,tr (Qh K h%)) : Q h 0,tr(Q) = l}. (14) 

32l as the joint field of values of the set of matrices h±hf , hxh^. The set 
1 The next result shows that the boundary of Q in any direction e £ £ can be 
achieved with rank-1 transmit covariance matrices. 
Lemma 3: Transmit covariance matrices from the set 

S := {Q : Q G 5, rank (Q) = l}, (15) 

achieve all points in d e Q, e G £. 

Proof: The proof is provided in Appendix ■ 
Lemma [3] supports the result in ll35l as a special case, where it is shown that single-stream beamforming 
is optimal in the MISO IC to achieve Pareto efficient operating points. Accordingly, efficient transmission 
strategies can be described by beamforming vectors. The next theorem characterizes the beamforming 
vectors that achieve the upper boundary of the set Q in a specific direction e. 

Theorem 1: All upper boundary points of the set SI in direction e G £ can be achieved by 

w(\) = v max h^hthf\ , (16) 



with 



AgA:=|ag [0,1]* :X> = lj. (17) 

Proof: The proof is provided in Appendix |D] ■ 
The interesting observation from Theorem Q] is that all upper boundary points in direction e of the K- 
dimensional gain-region can be achieved by a parametrization using K — 1 real parameters. Note that 
the components of the direction vector e in Theorem Q] can have positive and negative components. For 
the particular choices = for all I such that = +1, the largest eigenvalue of Z in ( fT6l ) is zero. 
In this case the largest eigenvalue of Z can also have geometric multiplicity larger than one, i.e., there 

The set is the convex hull of the joint numerical range of the matrices hih?,... , tiR-hx. The joint numerical range of 
a set of K matrices Ai, Ak is defined as 1331 : W(Ai, Ak) '■= {(z fl Aiz, z H Akz) : z H z = l}. In our case, 
Ai = h\hi , Ak = hxh^. Convexity of the set W(Ai, Ak) is not always satisfied. Conditions for its convexity are 
studied in 1341 . 
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exist multiple linearly independent eigenvectors associated with that eigenvalue. For example, assume 
N = 3, K = 2, e\ = +1, and e2 = —1- For Ai = and A2 = 1, the matrix Z in Theorem Q] is equal to 
—h,2h,2 which is rank-1 with the largest eigenvalue having geometric multiplicity two. In this case, we 
choose an eigenvector which lies in the span of H = [hi, Hk\. 

In Fig. |3l a three dimensional gain-region is plotted where the transmitter uses three transmit antennas. 
The boundary points are calculated by generating the beamforming vectors characterized in Theorem Q] 
For each boundary part, two real-valued parameters are required, and these are varied between zero and 
one with a step-length of 0.02. The gain-region is shown to have a convex boundary. Moreover, since it 
is possible to null out the power gain at two receivers simultaneously, the boundary of the gain-region 
touches each axis in one point. 

B. The number of transmit antennas satisfies N < K 

In case the number of antennas at the transmitter satisfies < ^T, it is not possible for the transmitter 
to choose a full power transmission strategy which nulls out the power gain at all except one receiver. 
Hence, the transmitter's ZF strategy would be to allocate no transmit power. This reveals that in order 
to achieve all boundary points of the power gain-region the transmitter has to reduce its transmission 
power. This is the reason why we study this case separately. 

We start by assuming that the transmit covariance matrices are chosen from the set S in ©. The 
corresponding power gain-region is compact and convex according to Lemma Q] Hence the following 
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programming problem, similar to the one formulated in (I35T ) in Appendix O achieves the boundary points 
of f2 in direction e £ £ : 

A" 

maximize X^e^hf Qh^ 

i=i (18) 

subject to tr (Q) < 1, Q y 0, 

where A € A is defined in ( fTTT ). As in Appendix [Cl the problem in ( fT8T ) can be equivalently written as 

AT / K \ 



maximize {j,g(Q)fj,t [ X^h^hj 



z (19) 

TV 

subject to < 1, ne(Q)>0, for all £ = 1, ...,iV 

£=i 

The solution of this problem is hn{Q) = P,P € [0)1]) ar *d w(Q) = for £ 7^ A. Hence, the 
optimal transmit covariance matrices Q are rank- 1 . In addition, the optimal power allocation is to allocate 
power only in the direction of the eigenvector corresponding to the largest eigenvalue of Z. Hence the 
formulation in Theorem Q] is still valid to determine the beamforming vectors that achieve the boundary 
points of fl However, power control is to be applied on specific beamforming vectors. In order to 
maximize the problem in (fl9l ), we choose p = 1 if the largest eigenvalue of Z is strictly larger than 
zero, and p = if the largest eigenvalue of Z is strictly less than zero. It can be easily checked that 
Z can be negative definite according to Weyl's eigenvalue inequality in ll36l Theorem 4.3.7]. We give 
a simple example to illustrate this case. Assume N = 2, K = 3, e\ = +1, and e2 = = —1. For 
Ai = and A2,A3 > 0, then Z = — A2ft.2^ — ^3^3^ which is negative definite. In case the largest 
eigenvalue of Z is zero, all feasible power allocations should be adopted in order to achieve the boundary 
of the gain region. For the same example given above and for fJ,(Z) = 0, the weighted sum gains in 
(|39l from Appendix |Pl formulates to X\X\{w\Wi) = X2X2{wiWi) + Xj,xj,{wiWi). This equation, with 
the variables in A, describes the plane that touches the boundary of the power gain-region at x(wiwf ) 
and passes through the origin. Choosing a power allocation p € [0, 1] achieves points on the segment 
connecting the origin and x(wiwf) which correspond to boundary points on the power gain-region. 
The power allocations that achieve the boundary points of the power gain-region U in direction e € £ 
are summarized as follows: 

p = 1 for fj, N (Z) > 0, p & [0, 1] for /J, N (Z) = 0, p = for fi N (Z) < 0. (20) 

where Z is given in ( fT9l ). 
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2 S 




Fig. 4. An illustration of a three dimensional gain-region where the transmitter uses two antennas. 

In Fig. HI a three dimensional gain-region is plotted where two antennas are utilized at the transmitter. In 
comparison to the plot in Fig. [3l the region looks like a cone whose vertex is at the origin. The outermost 
boundary which is not flat is attained with full power transmission. The flat surfaces correspond to the 
case for which the transmission power is varied between zero and one. 

IV. Pareto Boundary Characterization 

In this section, all T transmitters are considered again. The analysis of the single transmitter case in the 
previous section builds the framework to determine the beamforming vectors for each transmitter that are 
relevant for Pareto optimal operation. Each transmitter k, k G T, is associated with a power gain-region 
Slfc. The sets of intended and unintended receivers of a transmitter k are JC(k) and JC(k), respectively. In 
a general network, efficient operation of the systems requires the transmitters to maximize the power gain 
at intended receivers and simultaneously minimize the power gain at unintended receivers. In order to 
achieve Pareto optimal operating points a combination of these objectives is required at each transmitter. 
The gain-region characterization in the previous section, illustrates the effects of the beamforming vectors 
on the power gains achieved jointly at the receivers. Next, we formalize the transmit beamforming vectors 
of each transmitter that are relevant to achieve Pareto optimal points in the utility region. 
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Theorem 2: All Pareto optimal points in the utility region U can be achieved by beamforming vectors 

w k {Xk) = PkVmax ^2 ^k,eek,ih k: ih^ e , (21) 



with Afc G A defined in (fTTT ) and 



+1 t G K{k) 
-lie K{k) 



(22) 



i m(z k ) > o 

/'/, - S [0, 1] Miv(^fe) = • (23) 

o m(z k ) < o 

Proof: The proof is provided in Appendix [0 ■ 
Theorem |2] characterizes the transmission strategies that achieve Pareto optimal points. The number 
of real-valued parameters, A and p, that are required to characterize these beamforming vectors is: (i) 
T(K — 1) in case no power control is needed. This correspond to the case when N k > \JC(k)\. (ii) TK 
in case power control is needed corresponding to the case N k < \JC(k)\. All these parameters take values 
between zero and one. The direction vector in (l22l) can be determined since each transmitter k knows 
its intended and unintended receiver sets, K.(k) and K{k) respectively. This direction vector specifies the 
relevant boundary of the transmitter's power gain-region. 



A. Example Revisited 

The example in Section III-BI will be continued to find the beamforming vectors that achieve Pareto 
optimal points in the corresponding utility region. The choice of the beamforming vectors with unit norm 
for transmitters 11, 12 and 2 which achieve Pareto optimal points are characterized in Theorem |2] Since 
the number of antennas at each transmitter is strictly larger than the number of unintended receivers, all 
transmitters should transmit at full power according to d23l ). The power allocation p\\ and p\i on the 
beamforming vectors W\\ and W\i, respectively, is varied such p\\ + p\2 = 1. The power allocation 
parameters pn and pxi can be expressed as q and 1 — q, respectively, with q G [0, 1]. The characterization 
in Theorem [2] leads to the following nonnegative real- valued parameters: 

• For transmitter 11: Aii,i, An j2 , An,3, with An,! + An >2 + An,3 = 1. 

• For transmitter 12: Ai 2 ,i, Ai 2 ,2, Ai 2 ,3, with Ai 2j i + Ai 2j2 + Ai 2j3 = 1. 
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Fig. 5. Pareto boundary of the utility region of the setting described in Section HV-AI with SNR = 15 dB and N = 3. 



• For transmitters 11 and 12: q 6 [0, 1]. 

• For transmitter 2: A 2 ,i, A 2 ,2, A 2j3 , with A 2 ,i + A 2;2 + A 2j3 = 1. 

All seven required parameters are in the interval [0, 1], and the plot in Fig. [5] is obtained by varying these 
in a grid with 0.05 step-length. The points obtained include the points that lie on the Pareto boundary 
of the utility region which satisfy Definition Q] Points corresponding to weak Pareto optimality are not 
achieved by the parametrization in Theorem |2] and thus not included in the plot. Weak Pareto optimality 
is defined as in Definition Q] except that the corresponding inequality is not strict. Weak Pareto optimal 
points complete the shape of the utility region by orthogonally projecting each Pareto optimal point in 
Fig. [5] onto the coordinate surfaces. 

In the next section, we discuss two special applications of the developed framework. 

V. Applications 

A. Multiple-input Single-output Interference Channel 

The i^-users MISO IC consists of K transmitter-receiver pairs, where each receiver has an intended 
transmitter while all other transmitters induce interference on this receiver. We consider single-user 
decoding, i.e., interference is treated as additive noise at each receiver. For a given set of beamforming 
vectors {wi, wjc\, the following rate is achievable at receiver k, by using codebooks approaching 
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Fig. 6. Pareto boundary of a three user MISO IC rate-region, where N — 3. In the figure on the left, SNR = -10 dB, and in 
the figure on the right, SNR = 30 dB. 

Gaussian ones: 



This utility function satisfies properties A and B in Assumption Q] which leads to the following receiver 
sets for each transmitter k: K.(k) = {k}, and JC(k) = 1C\{k}. All points on the Pareto boundary of the 
achievable rate-region of the MISO IC can be reached by beamforming vectors as given in Theorem [2] In 
lfT3ll , a characterization of the beamforming vectors that achieve the Pareto boundary of the achievable rate- 
region is provided by a complex linear combination of the MRT and ZF strategies. The parametrization 
in Theorem [2] is real-valued with the same number of parameters, thus of lower dimension. 

The Pareto boundary of a three user MISO IC rate -region is plotted in Fig. [6] (left) for SNR = — 10 
dB and in Fig. [6] (right) for SNR = 30 dB. The generated points correspond to the beamforming vectors 
characterized in Theorem |2l The real-valued parameters are varied in a 0.05 step-length. Since we are 
only interested in revealing the Pareto boundary of the achievable rate region, we do the following. We 
randomly choose ten thousand generated points and remove all points that are dominated by these. In 
other words, for a randomly chosen rate tuple, all points corresponding to joint rates less than the chosen 
one are removed. In addition, the algorithm provided in 071 is applied. The algorithm reduces the number 
of plotted points in removing the points that are not visible from the viewed angle of the figure. This 
algorithm further reduces the complexity of rendering the generated points. 
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Next, we identify special operation points of each transmitter. 

1 ) Maximum ratio transmission: Maximum ratio transmission for transmitter k corresponds obviously 
to the parameters Xk.e = 1 for k = I and X^/ = for k ^ i. This transmission strategy is the unique 
Nash equilibrium (NE) strategy of each transmitter k, which is the outcome of a noncooperative game 
between the users. 

2) Zero-forcing transmission: Zero-forcing transmission is characterized with the following lemma: 
Lemma 4: Zero Forcing transmission for transmitter k corresponds to the parameters Xk,e = for 

k = £ and X k)i > for k ^ t. 

Proof: The proof is provided in Appendix |F] ■ 
The two-user MISO IC special case has an appealing form in terms of the parametrization of the Pareto 
boundary of the achievable rate -region. The efficient beamforming vectors are a linear combination of 
the MRT and ZF strategies lfl3l Corollary 2]. These two strategies have the interpretation of a transmitter 
being either selfish or altruistic (H, |[3T1 . The efficient beamforming vectors of transmitter k £ {1, 2} are 
given as lPT3l Corollary 2] 

wA) = jgriftzjj-r (25) 



A fc < RT + (l-A fc )™f 

where A& € [0, 1], w^ RT = tt^t, and iu| F = [T7ji M r~n > k ^ £■ Next, we prove that the parametrization 

|| K k £ || 

in (|25l ) has the same set of strategies as in Theorem [2] for K = 2. For this case, the eigenvalue equation 
for the hermitian matrix in Theorem Q] is written as (Xihuhn — (1 — ^h^h^^l = M" u, l- This 
equation can be equivalently formulated to Xi\\hii\\ 2 ^j^^Wi — (l — Xi)\\hi2\\ 2 jjr^ w i = ^W\. Adding 

1 1 'll 1 1 1 || »ll2 || 

(1 — Ai)||ft.i2|| 2 if i on both sides of the equation leads to 

Ai||h u || 2 n hll + (l-A 1 )||^ 12 || 2 n^ i2 )^ 1 = (^ + (l-A 1 )||/n 2 || 2 )t ( ; 1 . (26) 

The LHS of the equation states that the principal eigenvector is a linear combination of its orthogonal 
projection on h\\ and the orthogonal projection onto the orthogonal complement of h\2- Since the largest 
eigenvalue fi is larger or equal to zero, the weight in the RHS of d26l ) is always positive. Hence, the 
optimal set of beamforming vectors can be equivalently characterized by 



B. Noncooperative Underlay Cognitive Radio 

In an underlay cognitive radio scenario, secondary users can share the communication resources 
with primary users under the condition of not imposing quality of service (QoS) degradation to the 
primary systems. A limited QoS degradation to the primary users is described by interference temperature 
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constraints (ITC) [38]. When no interference on the primary users is allowed, the constraint is said to 
be a null-shaping constraint. Motivated by the concept of null-shaping constraints IT391 . the next result 
gives an alternative characterization of efficient transmission strategies. In order to be able to fulfill the 
null-shaping constraints, the number of applied antennas at the transmitter has to be greater than or equal 
to the number of primary receivers. We assume there exists virtual single-antenna primary receivers, and 
consider the efficient design of the null-shaping constraints. The following result has been presented in 
BUI without proof. 

Corollary 1: Assume N k > K and define the matrix 

Zk(^k) = z i(^k)i---,z\jc(k)\(^k),z Nk _^ k ^ +1 (X k ),...,z Nk ^i(X k ) , (27) 

where 

Zi(A fc ) = vi ^k,ee k ^h ki h^\ , (28) 

with Afc € A defined in (fTTT ) and e k ,e defined in (T22)) . All points on the Pareto boundary of the utility 

region IA can be reached by beamforming vectors 

II7 i\ \hhi 

w k (X k ) = , £eJC(k). (29) 



u z k (\ k ) h ke 

Proof: The proof is provided in Appendix [Gj ■ 
In Corollary [U the design of the null-shaping constraints is given in (T27T ). and the efficient transmission 
strategies are given in d29l . Here, K — 1 null-shaping constraints are to be applied on each transmitter, 
and the number of required real-valued parameters is the same as in Theorem [2] Hence, the complexity 
of parameterizing the efficient beamforming vectors is the similar in Corollary Q] and Theorem [2] 



The form of the transmission strategy in d29b has a relevant interpretation. Consider a MISO IC 

hkk 

setting as in the previous section. Rewriting d29l for transmitter k gives w k (\ k ) — — Zfc(Afc) 



n z k (x k) hkk\ 



This transmission strategy is MRT that satisfies null-shaping constraints given in Z k (X k ). The MRT 
strategy is the unique noncooperative strategy of transmitter k, and corresponds to the NE strategy 
in game theoretical terms [4-1 ] . Through the design of the null-shaping constraints in Corollary [T] all 
Pareto optimal points of the utility region are characterized by transmission strategies that are in NE and 
satisfy the characterized null-shaping constraints. The interesting observations are as follows. Null-shaping 
constraints are sufficient to characterize the Pareto boundary of the MISO IC rate region. Moreover, given 
the null-shaping constraints, the transmitters are required to be noncooperative in order to achieve efficient 
operating points 1401 . El . This result can be exploited by clever secondary user selection algorithms. 
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VI. Extension to Multiple Antennas at the Receivers 

We discuss the extension to multiple-antennas at the receivers and study which settings can be applied 
to our framework. Assume that the number of antennas used by transmitter k is and the number of 
receive antennas at receiver £ is Rf. The channel matrix from transmitter k to receiver I is denoted by 
H M E C RkXNk . The signal after receive filtering at receiver £ is 

T 

ye = ^2 z f H kew k s k + zfn e , (30) 

k=l 

where w k and zi are transmit and receive beamforming vectors, respectively. Considering a single 
transmitter k as in Section UTTl the power gain achieved at receiver t is x^{w k , zg) = \zf Hf-ew^l 2 . 
Clearly, if the receive beamforming vectors do not depend on the transmit beamforming vector Wk, then 
the framework in this paper can be applied to the setting. For example, the receive beamforming vectors Z£ 
can be fixed prior to transmission as in singular value decomposition (SVD) based receive beamforming. 
Let the SVDj of H^e be Uhu^HuVh ■ The SVD based beamforming vector z^ vd depends only on 
Hm and is the left singular vector in V h m which corresponds to the largest singular value in ^H ke - 
The framework in this paper can also be applied when the receive beamforming vectors depend only on 
intended transmitters' beamforming vectors. Assume the intended transmitter of receiver £ is k, then the 
maximum ration combining (MRC) beamforming vector of receiver I is z™ rc (wk) = HkeWk- In this case, 
the transmit beamforming vectors of the transmitter are not coupled through the receive beamforming 
vectors. Hence, the power gain region of a single transmitter can be used as a tool to characterize its 
efficient beamforming vectors. In case one receiver has multiple intended transmitters as receiver 2 in 
the example in Section ITl-B I Receiver 2 uses 2:2,12 = -H\i2)2' u 'i2 to successively decode the signal from 
transmitter 12 and uses 2:2,2 = H22W2 to successively decode the signal from transmitter 2. The utility 
of receiver 2 is log 2 (l + \w^ 2 Hi 2 H uw^ 2 + \w 2 H22H22W2 1 2 ), where the power gains from the 
transmitters are not coupled by the transmit beamforming vectors. 

The framework in this paper will not apply if the receive beamforming vectors depend on sev- 
eral transmit beamforming vectors. An example is the linear minimum mean square error (MMSE) 
receiver Il42l . Assume the intended transmitter of receiver £ is k, then the linear MMSE receiver is 
zf mse (wi, . . . , Wt) = ((J 2 1 + Ylj^k H jiw jwj H^)^ 1 H kiWk- Thus the transmit beamforming vectors 

5 H = Uh^hVh where Uh is the Rk x unitary matrix having the left singular vectors, Vh is the x unitary 
matrix having the right singular vectors, and Sh is the Rk x iVfc diagonal matrix containing the nonnegative singular values. 
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wi, . . . , wt are coupled through the receive beamforming vector. Hence, the power gain-regions of the 
transmitters are coupled by their beamforming vectors. 

VII. Conclusions 

In this work, we consider MISO wireless interference networks in which T transmitters and K receivers 
share the same spectral band. We provide a framework to determine the transmission strategies that 
are relevant for Pareto optimal operation. By studying the single transmitter's power gain-region, the 
properties of efficient transmission are acquired. We prove that the boundary of the gain-region is convex 
and achieved with single-stream beamforming. Due to the convexity of the boundary of the gain-region, 
the efficient transmit beamforming vectors can be parameterized by real-valued parameters. We determine 
and distinguish the conditions under which power control is required. When the number of antennas at 
the transmitter is greater than or equal to the number of existing receivers, we show that full power 
transmission achieves all boundary points of the gain region. In this case, the parameterizations of efficient 
beamforming vectors requires K — 1 real-valued parameters between zero and one. When the number of 
antennas at the transmitter is strictly less than the number of receivers, we characterize the transmission 
strategies for which power control is required. For this case, an additional real- valued parameter is needed 
that varies the power level at the transmitter. We apply the single-transmitter framework to the multiple- 
transmitter case. On determining the important boundary part of each transmitter's gain-region, all Pareto 
efficient beamforming vectors are characterized. This parameterizations simplifies the design of Pareto 
efficient resource allocation schemes. 
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Appendix A 
Proof of Lemma Q] 

Note that £1 does not correspond to the joint field of values |32] which is known to be convex and 
compact. Therefore, we provide a proof of the convexity of f2 and show that it is bounded and closed, 
thus compact. It is simple to show that Q is bounded because the power gain at the Ith receiver has a 
finite maximum which is achieved when the transmitter performs MRT to that receiver, i.e., x^(Q) < 
h e H £ - j = \\hi\\ . Therefore, the box described by the set y := <.x : < X£ < \\hi\\ j, contains Q, 
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i.e., U C y as illustrated in Fig. [2] The set Q is closed because the feasible set of transmission strategies, 
Q y and tr (Q) < 1, is compact and convex. Since every pre-image of a closed set is closed for 
continuous functions |43l , follows that O is a closed set. 

It remains to prove that is convex. For any two points x(Q x ) G Q and x(Q y ) 6 0, we prove that 
x(Q z ) G f2, where x(Q z ) = tx(Q x ) + (1 — t)x(Q y ) and i € [0, 1]. Any component of x(Q z ) is 

x e (Q z ) = tx e (Q x ) + (1 - t)x e (Q y ) = thfQ x h, + (1 - t)hfQ y h t = hf(tQ x + (1 - t)Qjfr*. (31) 

Hence, the transmit covariance matrices that achieve the line segment between x(Q x ) and x{Q y ) are 
given as Q z (t) = tQ x + (1 — t)Q y . Accordingly, since Q x and Q y are positive semidefinite, then Q z (t) 
is positive semidefinite for all t G [0, 1]. In addition, Q z {t) fulfills the trace constraint, tr (Q z (t)) < 1, for 
all t G [0, 1] since tr (Q z ) = tr (tQ x + (1 - t)Q y ) = ftr (Q x ) + (1 - i)tr (Q y ) < 1. Therefore, x(Q z (t)) 
also lies in Q, for all t £ [0, 1]. Hence, the set 17 is a compact and convex set. 

Appendix B 
Proof of Lemma [2] 

In order to prove that the transmit covariance matrices from S in (fb2l achieve points on the boundary 
of the set Q, in direction e, e E £, we show that for any transmit covariance matrix P with tr (P) < 1, a 
transmit covariance matrix Q from S can be constructed in which x(Q) dominates x(P) in direction e 
according to Definition [3] Assume ee = +1, we can construct Q as Q = P + (1 — tr (P))- ' 



\Xl^hthf Uz\\ ' 

where Z = [hi, /ty+i, fojf]. Clearly, Q is in the set S. Since K < N, the dimension of the 
null space of Z is greater or equal to one, therefore the projection Tig is not equal to the zero vector. 
The power gain achieved with Q at the fcth receiver, k G K,, is 

s*(Q) = fcf f P + (1 - tr (P)) "f^f"0 h*. (32) 

We distinguish two cases. For k ^ I, the power gain at the /cth receiver is 

x fc (Q) = hf (p + (1 - tr (P)) "^f"0 fc* 

y \\n^hih"n z \\ J 

njshthfni . (33) 



fcf Phfc + (1 - tr (P)) frf „f, £ ,i„ Z fefc = x k (P). 



n^hfUz 
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This implies that the gain at receiver k ^ I has not changed. For k = £, the power gain at receiver £ is 
Xe (Q) = hf ( P + (1 - tr (P)) "^"M hi = hfPhi + (1 - tr 



hf Ph, + (1 - tr (P)) ' I y\ > xe{p y 
\\U^h e h" U z \\ 



(34) 

According to the above results we can construct Q,Q G 5, such that x(Q) > e x(P) for any given 
P, with tr (P) < 1. Therefore, the boundary set <9 e O, e G f , can be achieved with transmit covariance 
matrices that fulfill the total power constraint with equality. 

Appendix C 
Proof of Lemma [3] 

Since the set Q is convex and compact, the boundary in direction e, <9 e S7, e G £ , can be achieved 
using the Supporting Hyperplane theorem |[29l Theorem 1.5] by the following programming problem 

K 

maximize Xie^hf Qhi 

i=i (35) 

subject to tr (Q) = 1, Q y 0, 

where A G A defined in ( fTTl ). The objective in 051 can be written as 

X if / K \ JV 

V XtethfQhz = V A^e £ tr (Qh*hf ) = tr Q V \te e h f hf = tr (QZ) < V fi t (Q)fi t (Z) 

ti ti \ tl J ti (36) 



z 

where the last inequality holds according to the von Neumann trace inequality of product of matrices 
1441 . Define Q = U qTiqUq and Z = U z^zU% where Uq and U z are unitary matrices. Then, 
the upper bound in (l36l ) can be achieved by choosing Uq = U z- Hence, the problem in (|35T ) can be 
equivalently written as 

N 

maximize ^ He(Q) He(Z) 
t=i 

(37) 

N v ' 

subject to V] MQ) = !> MQ) > 0, for all ^ = 1, iV. 

e=i 

The solution of this problem is /ijv(Q) = 1 and fie(Q) = if £ ^ N. Thus, Q is rank-1 and the transmit 
covariance matrices from S achieve the boundary of the region fl 
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Appendix D 
Proof of Theorem Q] 

As in the proof of Lemma [3j the convex boundary of the set 17 in direction e can be characterized by 
the solution of the following programming problem 



K 

2 



maximize X^et | w H h^ | 



1 1 2 

subject to ||tu|| = 1, 



(38) 



The objective function in (1381) can be rewritten as 

K / K \ 

y(w) = Y,^ e \w H h e \ 2 = w H lj2 X e e t h * h f ) w - ^ 39 ) 
l=x \i=\ ) 

S v ' 

z 

Note that the matrix Z in d39l is not necessarily positive semidefinite because the directional vector 
e can contain negative components. However, it is Hermitian and therefore, the solution to (l38l) is the 
eigenvector which corresponds to the largest eigenvalue of Z. 

Appendix E 
Proof of Theorem [2] 

The proof is by contradiction. Assume that a power gain tuple x k (w k w^) is not on the boundary of the 
power gain region £l k and achieves a Pareto optimal point in the utility region. Since Xj,(wj,w^) is not 
on the boundary of 17^, we can find another power gain tuple on the boundary of the power gain region 
which dominates Xk(wi,w^) in the direction specified in (1221) . Accordingly, it is possible to increase the 
utility of at least one receiver without changing the utilities of the other receivers which follows from 
Assumption [TJ Thus, a contradiction is made on the Pareto optimality of x^w^w^). Formally, assume 

x k {w k w% ) i d e *n k , (40) 

and achieves a point on the Pareto boundary of U, i.e. 

{uxixx^wxwf ), ...,x T ,i{w T w^)),...,u K {xi t K{wiWi), ■-, x t ,k(w t w% ))) G VB, (41) 

according to Definition [TJ By assuming d40b , we can find another feasible beamforming vector such 
that v k v^ € S k and achieves power gain tuple on the boundary of the power-gain region in direction 
efc such that 

x ki t(w k w%)e ki i < x k/ (v k v^)e k> i, for all £ e K,. (42) 
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Next we distinguish two cases for the inequality in (l42l) corresponding to I € K{k) or £ G /C(fc). 

1) Assume the inequality in (l42l is strict for £ G /C(fc) with = +1 as given in (1221 ). then 
Xk,i{ w k w k) < %k t( v k v k)- T ne P ower gains to all other receivers are to stay unchanged such that 

x k>j (w k w%) = x k ,j{v k v%) for j / £. Then, 

mix^eiwiwi), ...,x T ,e{w T WT)) < u^{xi >t {wxw" ), ...,x k> e(v k v k ), ...,xt,i(w t w^)) (43) 

holds according to property A in Assumption Q] in section III- A I This result contradicts (l4~Tl ). 

2) Assume the inequality in (l42l is strict for j G /C(&) with = —1. Then, (l42l changes to 

x k ,j{w k ) > x k j(v k ). Assuming x k j(w k ) = x k/ (v k ) for £ / j, then 

^•(xijftwiiwf ), ...,x T ,j(w T w^)) < u j (x 1 j(w 1 w^ r ),..., x kjj (v k v k ),..., x T ,j{w T wif)) (44) 

holds according to property B in Assumption Q] in section III-Ai This result contradicts (|4T1 ). 
Therefore, all points on the Pareto boundary of the utility region are achieved if each transmitter chooses 
its beamforming vectors to achieve the boundary of the gain-region in the direction specified in (1221 ). 
These beamforming vectors are characterized in Theorem [J in Section IIII-AI and for the optimal power 
control in (T20b in Section IIII-BI which lead to the formulation in Theorem [2] 



The zero-forcing condition on the beamforming vector of transmitter k, denoted as w£ , k G IC, is 



It is possible to fulfill the condition in (1431 ) if N k > K and all channel vectors from transmitter k to 
all receivers are linearly independent. Here, we give a direct proof. Assume the conditions X k ^ = for 
k = £ and > for k ^ £ hold, we have to show that these conditions lead to beamforming vectors 
that satisfy (l45l ). In order to do this, we have to study the matrix in (1211 ) whose eigenvector, corresponding 
to the largest eigenvalue, determines the used beamforming vector. Define the matrix M k , as 



B k 

where the direction vector e is specified as for the MISO IC application. According to the conditions 
X k ,£ = for k = £ and X k ^ > for k ^ £, A k in (1461 ) is equal to zero, and hence M k is 



Appendix F 



Proof of Lemma |4] 




(45) 




H 

■k(i 



(46) 
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negative semidefinite. The largest eigenvalue of Mk is therefore zero since Nk > K, and we can write 
( H ^k,ehkehke)iVk = 0- Since all channel vectors are linearly independent, the eigenvector associated 

tfk 

with the largest eigenvalue produces zero gain on any of the interference channel vectors. It is clearly 
seen that if \g = for any i ^ k, then the largest eigenvector does not necessarily produce zero gain on 
this receiver. 

Appendix G 
Proof of Corollary Q] 

We prove that the gains achieved by the beamforming vectors in ([291 are equal to the gains achived 
by the beamforming vectors given in (12T1 ). Define the matrix Mk, with the direction vector ek,e given 
in (l22l ). as 

K 

Mk = ^2 ^k,£ e k,lhk,lhk,l = ^k,eh k eh ke + 22 ~^k,thkthki 

£=i i€K{k) teK(k) ( - 47 ' ) 

> v ' V v ' 

A k B k 

The matrices Mk, Ak and Bk are Hermitian matrices of size iV& x N^. The eigenvalues of Mk are 
real and we always consider them ordered in nondecreasing order, i.e., ^\(Mk) < ^(Mk) < ... < 
/ijv fc (Mk). Ak consists of sum of positive semidefinite matrices. Hence, Ak ^ and rank (Ak) < \K.(k) | , 
i.e., 

< V\-j£(k)\+i( A k) ^ ••• ^ VN k (A k ), and m(A k ) = ... = V Nk -\Tc(k)\( A k) = °- ( 48 ) 

B k consists of the sum of the negative of positive semidefinite matrices. Hence, Bk ^ and rank (Bk) < 
|/C(/c)|, which leads to the following properties on the eigenvalues: 

fJ-i(Bk) < - < ii\jc{ k )\(B k ) < 0, and li\K(k)\+i( B k) = ■■■ = VN k ( B k) = 0- (49) 

Next, we study the eigenvalues of Mk = Ak + Bk- According to Weyl's inequality of the eigenvalues 
of the sum of Hermitian matrices [36, Theorem 4.3.7] the following properties are gained: 

V Nk -\Tc(k)\( M k) < V Nk ~\K(k)\( A k) + VN k (B k ) = (50) 
V\K(k)\+i(Mk) > Hi(A k ) + V\)c(k)\+i(Bk) = (51) 

The eigenvalues of Mk are ordered in nondecreasing order. Therefore, the following eigenvalues of Mk 
are always equal to zero: = ••• = ^N k -\Tc{k)\ (Mk) = 0. In addition, the smallest |/C(fc)| 

eigenvalues of Mk are nonpositive. 
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If the dimension of space is larger than the number of receivers, i.e., N k > \JC\, then there would be 

at least N k — \K(k)\ — |/C(fc)| eigenvalues of M k that are zero. For the eigenvectors corresponding to 

those eigenvalues, the eigenvalue equation is written as ( ^ ^k,ehkeh k £ + Yl —^k,fhkeh^ e )vi = 0, 

eelc(k) ' ieK(k) 
for all i = \K(k)\ + 1, ...,N k - \X(k)\. Then, for all £ e K, 



h,eh M hg) Vi = 0, for all i = \JC(k)\ + l,...,N k - \K(k) 



(52) 



The set of eigenvectors of M k , {vx, ...,v jv fc }, form an orfhonormal set, i.e. ||uj|| = 1 for all i = 
1, Nk and vfvj = for i ^ j. Therefore, we can write Yle=i v e v ^ = I> which gives 

N k -1 

v Nk (\)v§ k (X) = I - 5^ v e (\)vf(\) = I-G k (\)G%(\) = n£ kW , (53) 

i=i 

where G k (X) = [i>i(A), v^ k _i(X)] . Let the matrix Z k (X) consist of the eigenvectors of Gfc(A) 
excluding the eigenvectors that satisfy d52l) . i.e., 

Z k W= t»i(A),...,«| £ ( fc) |(A) ) t; JVfc _| )C(fc) | +1 (A),... > i>jvi b _i(A) , (54) 

then for any g £ C^* we can write 



■ I 



h 



1:1 



U Z k (X) h kt 



h 



k( 



\ v N k Wv% k (\)h ke \ 



\g H v Nk (\)\ 2 , (55) 



where i £ /C. Hence, the same power gains are achieved with 
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